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Alexander Merle@ 



Max-Planck-Institut fur Kernphysik 
Postfach 10 39 80, D-69029 Heidelberg, Germany 

In this paper, I give a pedagogical discussion of the GSI anomaly. Using two different formula- 
tions, namely the intuitive Quantum Field Theory language of the second quantized picture as well 
as the language of amplitudes, I clear up the analogies and differences between the GSI anomaly 
and other processes (the Double Slit experiment using photons, e + e~ — > u + fi~ scattering, and 
charged pion decay). In both formulations, the conclusion is reached that the decay rate measured 
at GSI cannot oscillate if only Standard Model physics is involved and the initial hydrogen-like 
ion is no coherent superposition of more than one state (in case there is no new, yet unknown, 
mechanism at work). Furthermore, a discussion of the Quantum Beat phenomenon will be given, 
which is often assumed to be able to cause the observed oscillations. This is, however, not possible 
for a splitting in the final state only. 
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INTRODUCTION 



In the last months, a measurement of the lifetime of 
several highly charged ions with respect to electron cap- 
ture (EC) decays at GSI Darmstadt [l| has caused a lot of 
discussion: Instead of seeing only the exponential decay 
law, a superimposed oscillation has been observed. The 
cause of this phenomenon, often referred to as Darmstadt 
oscillations or GSI anomaly, is not yet clear and a huge 
debate arose whether it could be related to neutrino mix- 
7 




not |fe1ll'[^,M"2l'lSsi3,M 26, ^(if. Alter- 
native attempts for an explanation involve spin-rotation 
coupli ng [29L [30l [3lT the interference of the final 
states [33[ , or hyperfine excitation [34[ • From the experi- 
mental side, two test-experiments (with, however, differ- 
ent systematics (35[) have been performed [36ll37l]. 
Several times in this discussion, the analogy of the GSI- 
experimcnt to the famous historical Double Slit experi 
ment using, e.g., photons [38| has been drawn [l,[li" 
which also led to lively debates at several meetings 
In this article, I show that the intuitive Quantum Field 
Theory (QFT) formulation of the problem always leads 
to the correct result. As there is still a lot of discus- 
sion in part of the community, it might be useful to give 
one more detailed explanation of the Quantum Mechanics 
(QM) involved. This can be done best by presenting easy 
and familiar examples that are not necessarily directly 
related to the GSI anomaly but do involve the same log- 
ical steps and are not under dispute. To do this, I start 
with the superposition principle and discuss the Double 
Slit experiment with photons, e + e~ — > scatter- 
ing, and the experiment performed at GSI. Afterwards, 
the language of amplitudes is used to further justify the 
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QFT-treatment by carefully considering several cases, 
where 7r + decay serves as additional example before the 
considerations are applied to the GSI-experiment, too. 
Furthermore, the so-called Quantum Beats (4fj| will be 
discussed, a well-known phenomenon that could indeed 
cause oscillatory decay rates. This has been argued to 
cause the observation at GSI in several places (see, e.g., 
Refs. 0, [12, [Hj])- It will, however, be shown that this 
cannot cause the observed behavior if a splitting is only 
present in the final state. 

In the course of the paper, we will see that all three lan- 
guages naturally lead to the same result, namely that 
a splitting in the final state cannot explain the GSI 
anomaly. Depending on the field of the reader, one or 
the other part may be clearer, but in the end it turns 
out that the intuitive QFT picture is correct and in per- 
fect agreement with the results obtained using probabil- 
ity amplitudes or the language of Quantum Beats, which 
are just different formulations of the same basic princi- 
ples. 



2. THE QUANTUM FIELD THEORY 
FORMULATION OF THE PROBLEM 

The starting point for the discussion is the superposi- 
tion principle in QM. One common formulation is [4lj | : 
"When a process can happen in alternative ways, we add 
the amplitudes for each possible way." The problem in 
the interpretation arises in the term "alternative ways" , 
because it is not a priori clear what the word "way" ac- 
tually means, as well as in the word "process" , which 
exhibits similar ambiguities. 

Let us use the following terminology: Process means a re- 
action with a well-defined initial and final state, whereas 
way is a particular intermediate state of a process. E.g., 
the scattering reaction e + e~ — > is one single pro- 

cess, no matter by which way (7-, Z -, or iJ°-exchange 
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Category 


Double Slit 


e + e — » fi + fj, 


GSI experiment 


1 


No slit-monitoring at all 


e + e~ -collider 


N/A 


2A 


Monitoring & read out 


N/A 


GSI-like experiment with more kinematical accuracy 


2B 


Monitoring without read out 


N/A 


Actual GSI-experiment 



Table I: The classification of the three examples. 



at tree-level in the Standard Model (SM) of elementary- 
particles) it is mediated. Z° — * v e v e and Z° —> v^v^ are, 
however, two distinct processes. 

Using this terminology, the superposition principle can 
be formulated in the following way: 

1. If different ways lead from the same initial to the 
same final state in one particular process, then one 
has to add the respective partial amplitudes to ob- 
tain the total amplitude. The absolute square of 
this total amplitude is then proportional to the 
probability of the process to happen [coherent sum- 
mation). 

2. If a reaction leads to physically distinct final states, 
then one has to add the probabilities for the differ- 
ent processes {incoherent summation). 

If a certain situation belongs to category 1, an interfer- 
ence pattern will be visible (or oscillations, in case the 
interfering terms have different phases as functions of 
time), while if it belongs to category 2, there will be no 
interference. The remaining question is at which point 
the measurement comes in. This can be trivially said for 
point 2: Either the experimental apparatus is sufficiently 
good to distinguish between the different final states (2A) 
- then no summation whatsoever is necessary simply be- 
cause one can divide the data set into two (or more) , one 
for each of the different final states. If this is not the case 
(2B), the experiment will be able to lead to either of the 
final states, but one would not know which one had been 
the actual result - then one would simply have to add 
the probabilities for the different final states to occur in 
order to obtain the total probability. 
What if we do such a measurement for category I? If we 
can indeed distinguish several ways that a process can 
happen, then this has to be done by some measurement. 
Since this measurement then has selected one particular 
way, we have actually transformed a situation belonging 
to category 1 into a situation of category 2. However, 
then there would be no terms to interfere with - the in- 
terference would have been "killed" . 
Let us now turn to Table [H that illustrates how our three 
examples fit into the categories 1, 2 A, and 2B. These 
three cases will be discussed one by one in the following. 

2.1. The Double Slit experiment with photons 

This is the "classical" situation of an experiment that 
reveals the nature of QM. It has first been performed by 



Thomas Young [38j | and has later been the major exam- 
ple to illustrate the laws of QM. Its basic procedure is the 
following: Light emitted coherently by some source (e.g. 
a laser) hits a wall with two slits, both with widths com- 
parable to the wavelength of the light. If it hits a screen 
behind the wall, one will observe an interference pattern, 
as characteristic for wave-like objects (category 1). There 
is, however, the interpretation of light as photons, i.e., 
quanta of a well-defined energy. Naturally, one could ask 
which path such a photon has taken, meaning through 
which of the two slits it has travelled. The amazing 
observation is that, as soon as one can resolve this by 
monitoring the slits accurately enough, the interference 
pattern will vanish, no matter if one actually reads out 
the information of the monitoring (2A), or not (2B). The 
reason is that, regardless of using the information or not, 
the measurement itself has disturbed the QM process in 
a way that the interference pattern is destroyed [42j . 
The key point is that one cannot even say that the photon 
takes only one way: In the QM-formulation, amplitudes 
are added (and not probabilities), and hence the photon 
does not take one way or the other (and we simply sum 
over the results) , but it rather has a total amplitude that 
includes a partial amplitude to take way I as well as an- 
other partial amplitude way 2. By taking the absolute 
square of this sum of amplitudes, interference terms ap- 
pear. 

A QFT-formulation involving elementary fields only 
would be much more complicated: One would sum over 
the amplitudes for the photon to interact with each elec- 
tron and each quark in the matter the slits are made 
of, after having propagated to this particular particle 
and before further propagating to a certain point on the 
screen. Of course, by using an effective formulation of the 
theory, one can find a much more economical description 
and the easiest one is to simply comprise all possible in- 
teractions into two amplitudes, one for going through the 
first and one for going through the second slit. 
Let us go back to this effective formulation: If there is 
monitoring, one actually "kills" one of these two ampli- 
tudes, the other one remains, and the interference is de- 
stroyed. Whenever there has been such a measurement, 
the interference will vanishe. As we will see, the ques- 
tion is if in a certain situation a measurement has been 
performed (or is implicitly included in the process con- 
sidered), no matter if the corresponding information is 
read out, or not. 
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2.2. e e — > /i /i scattering at a collider 

Let us now consider the scattering of e + e~ to form a 
pair of muons. This is, differently from the Double Slit 
experiment, a fundamental process where only a small 
number of elementary particles is involved. If one wants 
to calculate the scattering probability, the amplitude for 
the process is again decisive. In the SM, there are only 
three possibilities for this process to happen at tree-level 
and in all three of them the e + e~ pair annihilates to 
some intermediate (virtual) boson which in the end de- 
cays again, but this time into a pair. The inter- 
mediate particle can either be a photon, a Z°-boson, or 
a Higgs scalar, see Fig. [TJ 

Here, we have three different ways to form the process. 
The difference to the Double Slit experiment, however, 
is that these three ways cannot be separated easily. In a 
real collider-experiment we are not able to say that the 
reaction e + e~ — * has taken place by the exchange 

of, e.g., a photon only, but it will always be the sum of 
the three diagrams (and a lot more, in case we include 
higher orders). Hence, this process will always fall into 
category 1 and interference terms will appear. 
This is an easy and familiar example for the appearance 
of interference terms in a real experiment. In the next 
paragraph it will be shown exactly what is different in 
the case of the GSI-experiment. 

2.3. The GSI experiment 

The remaining question is what the situation looks 
like for the GSI-experiment. Even though the QFT- 
calculation of what happens is pretty straightforward, 
fitting everything in the language used above might be a 
bit more subtle. We will, however, see in Sec. 3.2 that the 
formulation in terms of amplitudes additionally justifies 
the result obtained here. 

Let us at first consider the Feynman diagram of the pro- 
cess involved in the GSI-experiment in Fig. [2] 17]: Here, 
in the absence of extreme kinematics (meaning that the 
Q-value of the reaction is large enough so that all neu- 
trino mass eigenstates can be emitted), the neutrino is 
produced as electron neutrino. What happens to this 
neutrino? Since it is not detected, it escapes to infinity in 
the view of QFT (in the picture of second quantization) . 
Physically, it loses its coherence after some propagation 
distance and travels as a unique mass eigenstate. 
The key point is the following: Since the neutrino will not 
interact before it loses its coherence, it must be asymp- 
totically a mass eigenstate. This can be shown easily: 
The coherence length of a (relativistic) neutrino is given 
by |H 

L coh = 2\f2a x ■ — — , (1) 
(Am 2 )© 

where a x is the size of the neutrino wave packet, p is 
the momentum of the mean value neutrino momentum 



in the limit m„ = 0, and (Am 2 ) = 7.67 • 1CT 5 eV 2 01 
is the solar neutrino mass square difference as known 
from neutrino oscillation experiments. The question is 
how to obtain an estimate for a x : If the nucleus was 
inside a lattice, one could estimate a width like the typ- 
ical interatomic distance, a x ~ lA, which would lead 
to L co h ~ 2 • 10 8 m. Of course, this precision cannot 
be reached in the GSI experiment. However, at least 
during the electron cooling [45j . the nucleus will be lo- 
calized to some precision. Since the velocity of the nu- 
cleus is known, this information could in principle be 
extrapolated for each run. A fair estimate would then be 
the average distance between two electrons in the cool- 
ing process, which is roughly give n by 1/ \Jn ~ 0.1mm, 
where n is the electron density [46[ . This leads to a more 
realistic coherence length of L C oh ~ 2 • 10 14 m. The pes- 
simistic case, where a x is taken to be the approximate 
diameter 108.36 ui/tt [47} of the Experimental Storage 
Ring (ESR) produces L co h ~ 6 • 10 19 m. The mean free 
path of a neutrino in our galaxy, however, is roughly 
1 • 10 m (for an assumed matter density in the Milky 
Way of 1 • 10 -23 g/cm 3 ), so the assumption that the neu- 
trino does not interact before losing its coherence is com- 
pletely safe. 

Even if we do not know in which of the three mass eigen- 
states the neutrino actually is, we know that it has to 
be in one of them. This knowledge is somehow obtained 
"a posteriori" , since the mass eigenstate only reveals its 
identity after some propagation. But, by conservation of 
energy and momentum, one could treat the process as if 
the kinematical selection had already been present at the 
production point of the neutrino. This "measurement" 
is enforced by the physical conservation laws. 
An analogous reasoning is given by Feynman and 
Hibbs [12], using the example of neutron scattering: Neu- 
trons prepared to have all spin up scatter on a crystal. 
If one of the scattered neutrons turns out to have spin 
down, one knows by angular momentum conservation 
that it must have been scattered by a certain nucleus. 
In principle, by noting down the spin state of every nu- 
cleus in the crystal before and after the measurement, 
one could find the corresponding scattering partner of the 
neutron without disturbing it. No matter if this would 
be difficult practically, by a physical conservation law 
one knows that a particular scattering must have been 
present, even if the corresponding nucleus is not "read 
out" . Accordingly, the corresponding interference van- 
ishes and the neutrons that have spin down after the 
scattering come out diffusely in all directions. 
This can also be formulated in the language of wave 
packets: We have complete 4-momentum conservation 
for each single component (which is a plane wave!) of 
the wave packets, but if we consider the whole wave 
packet, its central momentum does not have to be con- 
served [48|,|49|. However, all the different components can 
produce both possible neutrino mass eigenstates, but for 
a certain kinematical configuration of parent and daugh- 
ter components only one of the mass eigenstates will ac- 
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Figure 1: The diagrams contributing to e + e~ — > n~ in the SM. 
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Figure 2: The Feynman diagram for the GSI experiment. 



tually be produced. 

The rest is easy: If the GSI experiment had infinite kine- 
matical precision, one could read out which of the mass 
eigenstates has been produced and it would clearly fall 
into category 2A. Since, however, this information is not 
read out but could in principle have been obtained (e.g. 
by detecting the escaping neutrino), the GSI experiment 
falls into category 2B and one has to sum over proba- 
bilities. This logic works because we know that the neu- 
trino is, after some propagation, no superposition of mass 
eigenstates anymore, but just one particular eigenstate 
with a completely fixed mass. 

A viewpoint closer to the amplitude formulation would 
be: If the neutrino finally interacts, it has to "decide" 
which mass eigenstate it has, even if it was a superpo- 
sition of several mass eigenstates before. This is then 
equivalent to the image of having produced one particu- 
lar mass eigenstate from the beginning on. 



3. AMPLITUDES - PROBABLY THE EASIEST 
LANGUAGE TO USE 

In this section, I use time-dependent amplitudes for 
the different basis states to describe another example, 
namely charged pion decay, which I compare then to neu- 
trino oscillations (with referring to the actual situation 
in the GSI-experiment). The logical steps needed to un- 
derstand the familiar example of pion decay are exactly 
the same as the ones needed to understand what is going 
on at GSI. This description is clear enough to account 
for very different situations and allows for an easy and 
nearly intuitive understanding of the various cases. Fur- 
thermore, it yields an a posteriori justification of the view 
used in the preceding section. 



3.1. Charged pion decay 

It is well-known that a charged pion (e.g. 7r + ) can de- 
cay into either a positron in combination with an electron 
neutrino, or into the corresponding pair of /j.-like parti- 
cles. Let us consider the case of a pure (and normal- 
ized) initial state pion \tt + ). As this state evolves with 
time (and is not monitored), it will become a coherent 
superposition of the parent-state, as well as all possible 
daughter states: 

k+(t)) = A^{t)\-K + ) + A ll {t)\li + v IJt )+A e {t)\e + v e ), (2) 

where all time-dependence is inside the partial ampli- 
tudes Ai- Of course, this state has to be normalized 
correctly: 

\AAt)\ 2 + \A,(t)\ 2 + \A e (t)\ 2 = l, (3) 

with Ar(0) = 1 and A^O) = A e (0) = 0. One can 
understand Eq. ([2]) in the following way: The state at 
time t is a coherent superposition of the basis states 
{|tt + ), l^Vfj,}, \e + is e )} with time-dependent coefficients. 
Note that the basis states are orthogonal. The outcome 
of a certain measurement is some state |\I/): All that a 
detector does is projecting on just this state Of 
course, different detectors will in general be described 
by projections on different |\E')'s, which is a reflection of 
the influence of the process of measurement on the mea- 
surement itself. If one wants to know the probability for 
measuring that particular state, one has to calculate it 
according to the standard formula, 

P(*) = K*|7T+(t))| 3 . (4) 

The question is what |^) looks like. To make that clear, 
let us discuss several cases: 

• The (trivial) case is that there has been no detec- 
tion at all: Then we have gained no information. 
This means that the projected state is just the time- 
evolved state itself (we do not know anything ex- 
cept for the time passed since the experiment has 
started), and we get 

\{^ + (t))\ 2 = \{n + (t)U + (t))\ 2 =l. (5) 

This result is trivial, since the probability for any- 
thing to happen must be equal to 1. 
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• The next situation is when our experimental appa- 
ratus can give us the information that the pion has 
decayed, but we do not know the final state. Then, 
it can be either \fi + v^) or \e + v e ) and we remain 
with a superposition of these two states. The only 
information that we have gained is that the ampli- 
tude for the initial pion to be still there is now zero, 
A-x = in Eq. Then, the properly normalized 
state |*) is 

|^ = A ll {t)\ f j.+ u ll )+A e {t)\e+u e ) 

y/\Mt)\ 2 + \Mt)\ 2 

The absolute value square of the corresponding pro- 
jection is 

m* + m 2 = \Mt)\ 2 +\Mt)\ 2 , (7) 

and if there is any oscillatory phase in the ampli- 
tudes, Ak(t) = Ak(t)e lulkt , it will have no effect due 
to the absolute values. 

• What if we know that the initial pion is still 
present? This sets A^(t) = A e {t) = 0, and |*) 
is just A-* (t)]^) / x/\Att (t)\ 2 . The projection gives 

|(*|7T+(t))| 2 = |XW| 2 , (8) 

which again does not oscillate. 

• If one particular final state, let us say \e + v e ), is de- 
tected, then An(t) — Afj,(t) = and we get another 
term free of oscillations: 

\mn+(t))\ 2 =\Mt)\ 2 - (9) 



3.2. Neutrino oscillations and the GSI-experiment 

Let us now turn to neutrino oscillations. Here, as we 
will see, a state like |W) in Eq. (fTU|) can indeed be phys- 
ical in some situations. To draw a clean analogy to the 
experiment done at GSI, we consider a hydrogen-like ion 
as initial state \M) that can decay to the state \Dv e ) 
via electron capture. Since there was an electron in the 
initial state, we know that the amplitude for producing 
the mass eigenstate \v.i) is just U e i- If there is no relative 
phase between the two mass eigenstates, the neutrino 
produced in the decay is exactly the particle that we 
call electron neutrino. In any case, due to different kine- 
matics, the two mass eigenstates will in general develop 
different phases in the time-evolution. This means that, 
in spite of the mixing matrix elements U e i being time- 
independent, there will be a phase between the two neu- 
trino mass eigenstates. Completely analogous to Eq. ([2]), 
the time-evolution of the initial state will be given by: 

|M(i)) = AM{t)\M) + U el Ai{t)\Dv l ) + U <i 2A2{t)\Dv 2 ), 

(11) 

with \A M (t)\ 2 + \U el A!(t)\ 2 + \U e2 A 2 {t)\ 2 = 1 and 
Am(0) = 1- We can immediately look at different cases: 

• The parent ion is seen in the experiment: This 
kills all daughter amplitudes, Ai :2 (t) = 0. The 
only remaining amplitude is .Ajw(i), very similar to 
Eq. ([5]). With the proper normalization for |\&) one 
gets no oscillation again: 

K*|M(t))| 2 = |^ M (t)| 2 (12) 



The question remains when we do get oscillations at all. 
The answer is: It depends on what our detector mea- 
sures. If, e.g., the detector measures not exactly the state 
or \e + v e ), but instead some (hypothetical) super- 
position (e.g., some quantum number which is not yet 
known, under which neither /i + nor e + is an eigenstate, 
but some superposition of them), then one could measure 
the following (correctly normalized!) state: 



I*) 



V2 



(10) 



The squared overlap is 



1 



l(*k + (i))l 2 = 2 [K(*)l 2 + \A e (t)\ 2 + 2Vl{A;(t)Mt))] 

where the 25ft (.4* (t)A e (t)) -piece will, in general, lead to 
oscillatory terms. What has been done differently than 
before? This time, we have done more than simply killing 
one or more amplitudes in Eq. @, and this is the cause 
of oscillations: Whenever we are in a situation, in which 
the state playing the role of \^f) in Eq. (fT0|) is physical, 
the corresponding projection will yield oscillatory terms. 
As we will see in a moment, this is exactly what happens 
in neutrino oscillations. 



• The next case corresponds to the GSI-experiment: 
One sees only the decay, but cannot tell which of 
the two neutrino mass eigenstates has been pro- 
duced. This leads to Am(£) = and one has to 
perform a projection on the state 



I*) 



UeMtyDv!) + U e2 A 2 {t)\Du 2 ) 



y/\U el A!(t)\ 2 + \U e2 A 2 (t)\ 2 
Doing this with \M(t)) from Eq. dTTJ) yields 



(13) 



|<*|M(i)>| 2 = 



\U e Mt)\ 2 -l + \U e2 A 2 {t)\ 2 -l 



VWeMf}\ 2 

\U e2 A 2 (t)\ 2 , 



\U e2 A 2 (t)\< 



(14) 



which exhibits no oscillations, but is rather an inco- 
herent sum over probabilities. This result is the jus- 
tification of the intuitive treatment in Sec. 2.3: The 
elementary QM-discussion using probability ampli- 
tudes gives us just the correct prescription for how 
to sum up the amplitudes for the final states. 

• The GSI-experiment with infinite kinematical pre- 
cision: In this case, one could actually distinguish 
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the states \Dv\) and \Dv 2 ). If one knows that \Dv\) 
is produced (e.g., by having very precise informa- 
tion about the kinematics), one will again have no 
oscillation, 



K*|M(t))| 2 = |^(t)[/ el | 
just as in Eq. ([9]). 



(15) 



These are in principle all cases that can appear. One 
can, however, have a closer look at the realistic situation 
in the GSI-experiment. Let us re-consider Eq. (fTTjl : In 
reality, the parent ion will be described by a wave packet 
with a finite size or, equivalently, a finite spreading in 
momentum space, due to the Heisenberg uncertainty re- 
lation. If this wave-packet is broad enough that each 
component can equivalently decay into \Dv\) or \Dv 2 ), 
then both of the corresponding amplitudes will actually 
have the same phase {A\{t) = A 2 {ty), since they have 
the same energy, and one can write Eq. (|lip as 

\M{t)) = A M {t)\M)+A{t) [U el \Dvi) + U e2 \Dis 2 )] . (16) 

-> . ' 

=\Du e ) 

Since the knowledge of the momentum of the parent ion 
is not accurate enough at the GSI-experiment to make a 
distinction between both final states \Duk), this is a real- 
istic situation. Of course, this does not at all change the 
above argumentation, since the final state |\&) will expe- 
rience the same modification. The neutrino produced is 
an electron-neutrino, as to be expected. 
The question remains, why some authors come to the 
conclusion that there should be oscillations? The answer 
is simple: If the correspondence between time-evolved 
initial state and detected state is wrong, then oscilla- 
tions may appear. As example, we will consider the sit- 
uation that the kinematics of the parent and daughter 
are fixed so tightly, that indeed the production ampli- 
tudes for \Dv\) and \Dv 2 ) are not equal. This would 
correspond to an extremely narrow wave packet in mo- 
mentum space. Let us, e.g, have in mind the extreme case 
when by kinematics only the production of v\ is possible. 
This is no problem in principle and we would be used to 
it if neutrinos had higher masses, so that the Q-value of 
the capture was only sufficient to produce the lightest 
neutrino mass eigenstates. If only the disappearance of 
the parent is seen, the corresponding state |\E'), which is 
detected, is given by Eq. (fl3|) (with A 2 (t) = in the ex- 
treme case, but anyway with A\{t) ^ A 2 (t)). The corre- 
sponding neutrino is, however, no electron-neutrino any- 
more (which would be [/eil^i) + C/ e 2 1 ^2 } , with the same 
phase for both states)! Indeed this is no surprise at all, 
since the kinematics in the situation considered is so tight 
that it changes the neutrino state which is emitted. This 
is a clear consequence of quantum mechanics, since for 
obtaining the necessary pre-knowledge (namely the very 
accurate information about the kinematics), one has to 
do a measurement that is precise enough to have an im- 
pact on the QM state. 



If one consideres the state from Eq. ([13)) as being the 
one emitted but then projects onto an electron neutrino 
state, oscillations will appear: 

\{D,v e \M(t))\ 2 = \{U* el {D Vl \ + U* e2 (Dv 2 \) ■ 
■(A M (t)\M) + U el A x {t)\Dv x ) + U e2 A 2 {t)\Dv 2 ))\ 2 = 
= |_4x (i)| 2 + |.4 2 (t)| 2 + 2^{Ax{t)Ai{t)). (17) 

This is, however, wrong: One has not used all the infor- 
mation that could in principle have been obtained! But 
Nature does not care about if one uses information or not, 
so this treatment does simply not correspond to what 
has happened in the actual experiment. The oscillations, 
however, only arise due to the incorrect projection, and 
have no physical meaning. 

The remaining question to obtain a complete understand- 
ing of the situation is if the neutrino that is emitted in the 
GSI-experiment oscillates. The answer is yes, of course. 
But to see that, we will have to modify our formalism a 
bit. Knowing that an electron neutrino has been emitted 
corresponds to Am (t) — in Eq. (|16|). and the remaining 
(normalized) state is: 



I*) = 



At) 

\A(t)\ 



[U el \Dv x ) + U e2 \Dv 2 )]. 



(18) 



Re-phasing this state and measuring the time from t on 
gives as initial state: 



I*) = £7 e i|IM) + U e2 \Dv 2 ). 



(19) 



This is the state which will undergo some evolution in 
time according to 

|*(0) =AW)U*.\Dv 1 ) +A' 2 (t')U e2 \Dv 2 ), (20) 



With lAWWe 



\A' 2 (t')U e2 \ 2 = 1 and ^i(O) = 



A' 2 (0) = 1. If we ask what happens to this neutrino 
if it is detected after some macroscopic distance, it is 
necessary to take into account what has happend to the 
daughter nucleus that has been produced together with 
the neutrino, due to entanglement. The daughter nu- 
cleus, which is accurately described by a wave packet, 
is detected, but not with sufficient kincmatical accuracy 
to distinguish the different components \D) of the wave 
packet. The effect of such a non-measurement is studied 
most easily in the density matrix formalism. The density 
matrix p' corresponding to Eq. (|20p is given by 



\^{t'))^{t')\ = \B 1 {t')\ 2 \D)\ Vl )( Vl \(D\ 
+ \B 2 (t')\ 2 \D)\v2)(v2\(D\ + 
+[B l {t')B* 2 tf)\D)\v l ){v 2 \{D\+h.c], 



(21) 



where Bt{t') — A' k (t')U e k- If the exact kinematics of the 
daughter is not measured, then one has to calculate the 
trace over the corresponding states. It gives 

p = J dD(D\p'\D) = |B 1 (i')|>i)(^i| + (22) 
+ \B 2 {t')\ 2 \v 2 ){u 2 \ + (B 1 (t')B* 2 (t')\iy 1 )(iy 2 \ + h.c). 
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Figure 3: Type I (left) and type II (right) of the Quantum 
Beats settings. 



If we want to know the probability to detect, e.g., a muon 
neutrino, \v^) = U^Vi) + U^ 2 \v 2 ), the corresponding 
projection operator is given by 

^ = K)(^|, (23) 

and the probability to detect this state is 

P M = Tr(7» = H7>h) + H7>h>. (24) 

Note, however, that the neutrino states \v\. 2 ) will always 
be orthogonal, since they correspond to eigenstates of dif- 
ferent masses (like an electron is in that sense orthogonal 
to a muon). The result is 

P, = \U, 1 \ 2 \B 1 (t')\ 2 + \U» 2 \ 2 \B 2 (t>)\ 2 + 

+ [U^U; 2 Bl(t')B 2 (t') + cc], (25) 

whose second line contains oscillatory contributions. 
These oscillation are indeed physical: Eq. (|23[) is a de- 
scription of a detector that is sensitive to z^'s only. If 
it could not distinguish different neutrino flavours, the 
oscillation would vanish again. 



4. QUANTUM BEATS 

The last point to discuss are the so-called Quantum 
Beats (QBs) [4(| ■ This phenomenon is known from Quan- 
tum Optics and has often been mentioned as possible 
explanation for the GSI anomaly. As we will see, the 
corresponding language can be equally used to describe 
the GSI-experiment and (of course) yields the same re- 
sult as already obtained. Still, it is also useful to consider 
the experiment from this point of view in order not to be 
misled by claims that erroneously make QBs arising from 
a splitting in the final state responsible for the observa- 
tion at GSI. 

Normally, one considers atomic levels for this discussion, 
and we will stick to that here for illustrative purposes and 
give the relation to the GSI-experiment at the end of each 
section. This way also easily clarifies the analogies to the 
Quantum Optics formulation. 



4.1. Single atom of type I 

Let us start with the classic example of QBs, namely an 
atom in a coherent superposition of three states \a), \b), 
and |c), where the first two states are above and closely 
spaced compared to |c). This setting is drawn on the left 
panel of Fig. [3] and is referred to as "type I" . First note 
that the three levels correspond to different (but fixed) 
eigenvalues of the energy and are hence orthogonal vec- 
tors in Hilbert space. This is not at all changed by an 
energy uncertainty which, however, makes it possible to 
have a coherent superposition of the three states. Ini- 
tially, we assume the atom to be in such a superposition 
of these states, but having emitted no photon yet. Ac- 
cordingly, the photon state can only be the vacuum |0) 7 . 
Then, the initial state of this system can be written as 



|*(0)) = A|a)|0) 7 + B„|b)|0) 7 +e o \c)\0) 



7! 



(26) 



where |.Ao| 2 + \Bo\ 2 + |C | 2 = 1. If this system undergoes 
a time-evolution, the lower state might be populated by 
de-excitation of the upper ones, which is done by photon 
emission. If the state |l x )7 = a £|0) 7 is assumed to de- 
scribe a state with one photon of frequency lo x , then the 
state at time t can be written as 

|*(t)) = A(t)\a)\0^ + B(t)|6>|0) 7 +C{t)\c)\Q) 1 + 
+C 1 {t)\c)\l ac ) 1 +C 2 {i)\c)\l bc ) 1 , (27) 

where A(0) = A , 8(0) = B , C(0) = Co, C 12 {0) = 0, 
and |„4(i)| 2 + \B(t)\ 2 + \C(t)\ 2 + |d(t)| 2 + \C 2 (t)\ 2 = 1. 
Under the assumption that all levels are equally pop- 
ulated, the radiated intensity will be proportional to 
(#(i)|E 2 (0,t)|#(i)), where 



E(x, t) = e k ,A (aic. 



-ikx 



at e+«* 



k.A 



(28) 



is the electric field operator and ek,A is the electric field 
per photon of momentum k and polarization A. Note that 
the creation and annihilation operators have only one 
non-trivial commutation relation, namely [at, a, cl^' a'] = 
<5k,k"5A,A'- In our case we obtain effectively: 

E(0,i) 2 = e 2 ac (l + 2ala ac )+ (29) 
+e 2 bc (l + 2al c a bc ) + 2e ac e bc (ai c a bc e iAt + a f bc a ac e- iAt ), 

where A = uj ac — Lj bc . Here, we have already used that 
terms like, e.g., a 2 c give no contribution with |\E') from 
Eq. I[27p. Remember now, that the atomic states are or- 
thonormal. This means that one can, e.g., combine a 
term proportional to (6| in ('F(i)l only with the corre- 
sponding term \b) in |\&(t)). The corresponding combi- 
nation of amplitudes |S(£)| 2 does, however, not oscillate, 
since any phase will be killed by the absolute value. This 
is also true for every term involving one of the constant 
parts of Eq. (|2"9"|) : E.g. the term proportional to C*{t)C\{t) 
can involve a factor 



|aLowaL|0> 7 = 0, 



(30) 
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because of a\ c acting on the left. There are, however, 
remaining oscillatory terms such as Cl(t)C2{t)e tAt , which 
is proportional to 



\a a , c a\ c a bc a\ c \[)) 1 = 7 (0|(l+4 c a ac )(l+4 c a fec )|0). 



1. 



These terms cause the Quantum Beats for a type I atom. 
Actually, one could have expected this result intuitively: 
Both of the coherently excited upper levels can decay into 
the same state |c) via the emission of a photon. Hence, 
one cannot in any way determine the photon energy with- 
out measuring it directly. Without such a measurement, 
interference terms will appear. 

How is the situation for the GSI-experiment? In this 
case one simply has to replace the photon by the neu- 
trino. As explained in Ref. [TtJ for instance, a split- 
ting in the initial state could lead to an oscillatory be- 
havior. This splitting, however, would have to be tiny, 
r~> 1CP 15 eV, a value which can hardly be explained. Fur- 
thermore, there exists preliminary data on the lifetimes 
of 142 Pm 60+ with resp ect to /3 + -decay that shows no os- 
cillatory behavior [12j | . An initial splitting in the nucleus 
would lead to an oscillatory rate in this case, too. Ac- 
cordingly, if such a splitting is present in the initial state, 
it could be in the levels of the single bound electron, since 
this would then affect EC-decays while leaving /3 + -decays 
untouched. 



4.2. Single atom of type II 

We can study a similar setting, namely an atom of 
type II, shown on the right panel of Fig. [3] The cor- 
responding initial state would again be described by 
Eq. (J2U), but its time-evolution would now look like 



mt))=A(t)\a)\0) 1+ B(t)\b)\0) 7 
+Z?'(t)|6)|l Qb ) 7 +C'(t)|c)|l QC ) 7 , 



■C(t)\c)\0) 1 + 
(31) 



where A(0) = A , 6(0) = B , C(0) = C , E'(0) = 0. 
C'(0) = 0, and \A(t)\ 2 + \B{t)\ 2 + \C{t)\ 2 + \B'l(t)\ 2 + 
|C'(i)| 2 = 1. The square of the electric field has again 
the form of Eq. (f2"5)) . just with be — > ab. Due to the 
orthogonality of the atomic states, there are not too many 
combinations which are possible: 

• 0-photon state coupled with itself: 
If we take, e.g., the term |„4(i)| 2 , it does not oscil- 
late anyway. Hence, only the time-dependent parts 
in Eq. ([25]) (with be — ► ab) could lead to oscillations. 
But they are proportional to 



|aLaah|0) 7 = 7 (0|<4 6 a ac |0) 



0. 



1-photon state coupled with itself: 
|£>'(t)| 2 does not oscillate, too, and the time- 
dependent terms from the electric field yield 

7 (l a6 |(4 c a Qf) |l a6 ) 7 = 1 {Q\a a ba\ c a ab a\ b \Q) 1 = and 
1 {l a b\al b a ac \l ab ) 1 = 7 (0| a a6 a^a ac o^ |0) 7 = 0, 



which follows immediately from the action of a\ c 
to the left and of a ac to the right, respectively. 

• 0-photon state coupled with 1-photon state: 
This is the only possibility, which is left. If we take 
for instance the term B*(t)B'(t), this will oscillate 
anyway, so we will also have to check the constant 
terms in Eq. (|2U)) . The ones proportional to 1 are 
naturally zero, 7 (0|l a &) 7 = 7(0|a^ b |0) 7 = 0. The 
other terms are 

7 (0| a\ c a ac |l a6 ) 7 = 0, 7 (0| a[ b a afc |l a6 ) 7 = 0, (32) 
7 (0| a) ac a afc |l ab ) 7 = 0, and 7 (0| a\ b a ac |l ah ) 7 = 0, 



where the action of the operators to give zero is 
always indicated by the arrow. The argumentation 
is analogous for the complex conjugated term. 

Hence, there can be no Quantum Beats for a single atom 
of type II! The intuitive reason is that, by waiting long 
enough, one could reach an accuracy in energy that is 
good enough to distinguish the possible final states \b) 
and |c). This would then be a way to determine the 
energy of the emitted photon without disturbing it. 

To give an analogous argument for the GSI- 
experiment, one has to turn the comparison given in 
Sec. 4.1 round and replace the atom by the neutrino and 
the photon by the ion. The reason is that what is claimed 
to interfere in this situation is the neutrino states them- 
selves (see, e.g., Ref. Q). This neutrino is not expected 
to interact, before losing its coherence (cf. Sec. 2.3). How- 
ever, once it interacts, it has to decide for a certain mass 
eigenstate. By monitoring this interaction, it would in 
principle be no problem to determine the neutrino's mass 
(e.g., by exploiting the spatial separation of the mass 
eigenstates far away from the source) and from this one 
could easily reconstruct the kinematics of the daughter 
ion in the GSI-experiment. Accordingly, no QBs are to 
be expected in this situation. 



4.3. Two atoms of type II 

On the other hand, there is a situation in which we 
can expect QBs even for atoms of type II, namely if we 
have two of them. If these two atoms are separated by 
a distance which is smaller than the wavelength of the 
emitted photons, there is no way to resolve their sepa- 
ration in space and we have to write down a combined 
initial state for both atoms, 1 and 2: 

|*(0)) = A|a)i|a) 2 |0) 7 +B |6)i|6) 2 |0) 7 +C |c)i|c) 2 |0) 7 
+^i,o|a)i|6)2|0) 7 +P 2 ,o|fe)i|a)2|0) 7 +£i,o|a)i|c) 2 |0) 7 + 
+f2,o|c)i|a) 2 |0) 7 + ^i,o|6)i|c) 2 |0) 7 +^2,o|c)i|6) 2 |0) 7 . 
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The corresponding time-evolution |\&(f)) looks a bit com- 
plicated: 

^(t)|a) 1 |o) 2 |0) 7 +B(t)|6)i|6)a|0) 7 +C(t)|c)i|c) 2 |0> 7 + 
+X>i(t)|o)i|6) 2 |0) 7 +Z>2(t)|6)i|a)a|0) 7 + 
+£i(t)|a)i|c) 2 |0) 7 + £ 2 (i)|c) 1 |a) 2 |0) 7 + 
+J 7 i(t)|6)i|c) 2 |0) 7 + ^ a (t)|c)i|6) 2 |0) 7 + 
+Si(*)|6)i|a) 2 |l at ) 7 + e 2 (t)|o)i|6) 2 |l a 6> 7 + 
+Wi(t)|c) 1 |a) 2 |l ac ) 7 +W 2 (t)|a) 1 |c) 2 |l ac ) 7 + 
+2k(i)|6)i|6)a|la6) 7 +Z 2 (t)|c) 1 |c) 2 |l oc ) 7 + 

+Jl(t)\b) 1 \c) 2 \lab)~ t + J2{t)\c) l \b) 2 \l ab ) 1 + 

+/Ci(t)|6)i|c> 2 |l„ c > 7 +/C 2 (i)|c)i|6) 2 |l ac ) 7 . (33) 

One oscillatory term would then be, e.g., J*K.\e~ %At , 
which is proportional to 

■y(^ab\a f ab a ac \l ac ) y = 1 {Q\a ab a\ b a ac a\ c \Q) 1 = (34) 

= 7 (0|(1 + a[ b a ab )(l + a{ c aac)\0} 7 = 7 (0|0) 7 = 1. 

0<- ->o 

If the spatial separation is less than the photon wave- 
length, one cannot determine the photon energy, because 
one docs not know which atom has emitted the radiation. 
Accordingly, we expect QBs. 

For the GSI-case, this possibility has to be taken into 
account, because even for runs with one single EC only, 
there can have been more ions in the ring that were lost 
or decayed via /3 + . In this case (comparing the neutrino 
again with the photon), one has to replace the wave- 
length of the photon by the de Broglie wavelength of the 
neutrino. The neutrino energy should be of the same or- 
der as the Q-value of the EC-reaction, which is roughly 
1 MeV [l]. The corresponding wavelength is, however, 



A = ~ 10~ 12 m, while the average distance between 
two ions should be of the order of the storage ring [5fJ , 
which is roughly 100 m [471 ] . Hence, this possibility is 
excluded for the GSI-experiment. 



5. CONCLUSIONS 

A comparison of the GSI-experiment with several other 
processes (the Double Slit experiment with photons, 
e+e~ — > [i + /x~ scattering, and charged pion decay) has 
been given. By using the language of QFT as well as 
the intuitive formulation with probability amplitudes, I 
have shown that the situation at GSI cannot lead to any 
oscillation of the decay rate, if the correct treatment is 
chosen and no additional assumptions (as, e.g., a split- 
ting in the initial state) are taken into account. Also 
the frequently mentioned possibility of Quantum Beats 
of the final state cannot explain the observed oscillations, 
at least not in the standard picture. Hopefully this ar- 
ticle will contribute to the clarification of the physical 
situation in the experiment that has been performed at 
GSI. 
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